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Introduction 

The formation of multiple wave fronts is important in applications of singularly 
perturbed parabolic systems. These solutions can be effectively constructed by 
formal asymptotic methods. When truncated to a certain order in e, they become 
formal approximations of solutions to the given system. The precision of a formal 
approximation is judged by the smallness of the residual error in each regular and 
singular layer and the jump error between adjacent layers. 

The purpose of this paper is to introduce Spatial Shadowing Lemmas that 
help to construct exact solutions near the formal approximations. The Shadow- 
ing Lemma was first developed for discrete mappings in K", see [8]. It has been 
extended to continuous flows governed by ODEs [11], and semiflows governed by 
abstract parabolic equations [2]. We will call these Temporal Shadowing Lemmas 
since the dynamical systems considered there evolve in time. 

The Temporal Shadowing Lemma has been used to construct exact solutions 
for singularly perturbed ODEs [11]. However, it cannot be applied directly to 
singularly perturbed parabolic equations. For formal approximations of multiple 
waves, the jumps between adjacent layers are functions of t and they occur along 
the x-direction. Since parabolic equations cannot be solved in the x-direction, 
therefore, they do not define a dynamical system in the spatial direction. 

To solve this problem, we use an idea motivated by the works of Kirchassner 
and Renardy [10, 15]. We find stable and unstable subspaces of the trace space 
so that the parabolic system can be solved forward and backward in the spatial 
direction. Thus the jumps along the lateral common boundaries can be corrected 
using the technique of the usual Shadowing Lemma in abstract spaces. 

This paper is divided into two parts. In the Part I, we show that for a general 
parabolic system, if a formal approximation is precise enough, then there is an 
exact solution near the formal approximation for at least a short time. The result 
obtained here applies to various systems including reaction-diffusion equations, 
Cahn-Hilliard equations [1], and viscous profile of conservation laws. In Part II, 
we show that with additional restrictions, the process in Part I can be repeated 
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to obtain global solutions if the formal approximation is a global one. Examples 
include reaction-diffusion equations and phase field equations [7] . 

PART I. Local Existence of Multiple Waves 

1. Consider a general singularly perturbed parabolic system, 

(1) eu t + (-e 2 ) m D 2 x m u = f(u,eu x ,--- , (eL^) 2 ™-^, x, e), u e R", x E K, 

where / is C°° with bounded derivatives in all the variables. Assume that the 
system has regular and internal layers located alternatively along the x axis. For 
simplicity, we solve the system for with no boundary conditions other than 

u e if 2 " 1 ' 1 . Assume that there are curves T' = {(x, t) : x = x l {t), t e [0, At}}, i e 
Z, that divide the domain x £ R, t £ [0, At] into regular or singular layers each 
is between T 1 ^ 1 and T l . 

Assume that a formal approximation is given, piecewise continuous, with 
u = w % {x, t, e) in S J . Using the stretched variables £ — x/e, t = t/e, let w J (£, t, e) = 
w l (e£, ct, e) which depends slowly on r. Assume that u> 1 € if 2m,1 (S J ). Let W l — 
(vf, D^w 1 , - ■ ■ , D 2m ~ 1 ) T . The error terms in the followings are — g l and —S l , 

(2) < + (-ir£>f"V - /(«;*, ■ ■ ■ , e^, e) = in £*, 

(3) TU J (f,r)-TU J+1 (f,r) = -^(T), at P. 

Let an exact solution be u l + w l in £\ Let U* = (u\ D 6 u l , • • • , D 2 " 1 " V) T . 
Let r e I = [0, Ar] where At is independent of e. The domain [0, At] corresponds 
to a short time interval [0, eAT] in the t variable. Assume that e is small so that a 
near identity change of coordinates in Y? can be used to straighten the boundaries 
T 1 ^ 1 and T l . In the following, we assume that £ l = x l /e is independent of t, but 
may depend on e, ft 1 — and £* = O* x i. With the new coordinates 

introduced above, linearizing (1) around at the fixed time r = 0, we have 

2m-l 

(4) < + (-lr^DfV - ^ 4(0^ =M\u\ g\ e), in E*, 

(5) [/ J (r,r)-C/ J + 1 (r,r)= ( 5 J (T), at P, 

(6) u%0) = ui(0, mSl\ 

Here ^4}(£) = Djf(w l (£, 0, e), D^w l (^, 0, e), • • • ) is the partial derivative of / with 
respect to the j-th variable. u l E H m (Q l ). N l depends slowly on t and I^Il 2 ^) 
= 0(\u% 2m , 1{j:i) + \g%, + leATlluV-^E'))- 

We look for a sequence of solutions {V}^ with u l E if 2 " 1 ' 1 ^). Let H k (I) 
be the usual Sobolev space and let Hq(I) be the completion of C°° functions, 
which are zero in a neighborhood of t — 0, in the H k norm. Define the product 
spaces 

B m (I) = n 2 ™-^ 1 "^^), B^{I) = Ul^H^il). 
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From the Trace Theorem, [14], the mapping 

£->tf(£,-), W^B m (I), 

is continuous with 

\U\t-)\B^ { i)<CW\ H , m ^^ 
Therefore, 5 l e B m {I) in (3) and (5). Let U l = (4, • • • 



,£)™" 1 4) r . Let TTj be the 

projection to the first j-tuples in the product space IL^R™, i.e., ttj (ui , • • • , U2m) = 



(til, ■ ■ ■ , tij), A compatibility condition is also assumed on the initial data and the 
jumps, 



(7) 



L73 +1 (f 



»**(<>). 



Notice that only the first m components of 5 l have well defined traces at r = 0. 

We now discuss the method of solving (4)-(6) with the compatibility (7) in 
Sections 2-5. 



2. Consider 

with it(£,0) 
system, 

(8) 



u T + {-l) m Dl m u = 0, £eR, rel+, 
0. Applying the Laplace transform, we have the so called dual 



U t = J(s)U 



I 
I 



\(-l) m+1 s 





0/ 



u. 



The matrix J(s) has 2m eigenvalues A = [(— l) m+1 s] i . 
Consider a sector in C, 

5 e (M) = {s : |s| > M, \arg(s)\ <6}, M > 0, ir/2<9 < n. 

When s <E Sg(M), the eigenvalues A are in 2m disjoint sectors of C, with |RcA| > 
cos (lim^ ) 2i \/N"- There are m eigenvalues with positive real parts and m with neg- 
ative real parts. Each eigenvalue has an n-dimcnsional eigenspace spanned by 
(u, Xu, • • • , \ 2m ~ 1 u) T , u e K". 

Let E m ' v (s) be the Banach space of points in R 2nm with an s-dependent 
norm, 



(U ,Ui, ■ ■ ■ , U2 m -l)\E™."(s) 



2m -1 
3=0 



|M2m-l-j| 



We actually will only use v = or j!- in this paper. Let P s and P„ be the 



projections in 



»2mn 



to the stable and unstable spaces of J(s), s G Sg{M). The 



projections are of rank ran and can be constructed using eigenvalues and eigen- 
vectors. Using the E m ' u {s) norm, we can show that there exist K, a\, a > 0, such 
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that 

\e J ^P s \ Em ,, {s) <if e - ai2 VN« <ic e -tt(i+ 2m /NK ] ^> 0j 
|e J W«P„| sm „ (s) <Xe« l2 V^« < ATe«( 1 + 2 \^)« , £ < 0. 
3. Consider 

2m-l 

(10) UT + (-ir£>| m «- ^ A)(ODiu = o, teg- 1 ,?], 

3 = 1 

with 0) = 0. Using the Laplace transform, we have a dual system, 

/ \ 

(11) U t = J(a)U+(-ir\ ••• ]U. 

Let T l (£> £, s) be the solution matrix for system (11). System (11) is said to have 
an exponential dichotomy in E m ' v (s) for £ e [£ l_1 >£ l ]i s € S$(M) if there exist 
projections Pj(£, s) + P„(£, s) = 7 in E m '"(s), continuous in £ and analytic in s, 
and positive constants K, a, such that 

|r(£,C,s)Pi(C,s)| E —( s) <Ke-^ 1+2m V^^\, £>C, 

|r(C,C,^(C^)|^ (s) <if e -«d+ 2 VN)i?-ci e < c . 

Suppose that sup ?G [ a h ] 0<fc<2m _i 1^(01 — From the Roughness of Expo- 
nential Dichotomy Theorem, [3], which is also valid in the Banach space E m - V (s) 1 
we find that there exists M = M(C) > 0, sufficiently large, such that (11) has an 
exponential dichotomy in E m - V (s) for £ e [f ~\f ] and s <E 5 e (M). On the other 
hand, if M is fixed, then there exists C = C(M) > 0, sufficiently small, such that 
the same conclusion holds. 

A function f(s) is in the Hardy-Lebesgue class ^(7), 7 € K, if 

(i) f(s) is analytic in Re(s) > 7; 

(ii) {sup^a^ \f (a + iu,)\*du>)Y'* < co. 

H(j) is a Banach space with the norm defined by the left side of (ii). Based on the 
Paley- Wiener Theorem, [16], if e~^f{t) e 7 2 (M+), then f(s) G ^(7), vice versa. 
For ft > and 76ft, define a Banach space 

H k (j) = {u{s) I u(s) and (s - j) k u(s) e W( 7 )}, 

M« fc ( 7 ) = Mwb) + l( s - 7) fc w|w( 7 )- 
For any 7 e R, ft > 0, there exists C = (7(7, ft) such that 

C- 1 (l + \a\ k ) < l + |s-7| fe < C(l + |s| fe ). 
Therefore an equivalent norm for H k ('j) is 

\u(<7 + iuj)\ 2 (l + |cr + iuj\ 2k )duj. 

-00 

It can be shown that if e"T*/(i) e H k {R+), then /(s) € ^(7). 
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We often use Banach spaces of functions with norms weighted by e 7T . For 
example B% l (R+,~f) = {4> : e^ 7 ~4>(t) e Bo"(R+)} with obvious norms. Other 
spaces like L 2 (R x R + ,7), etc. can be defined similarly. 

Let /C m (7) = n^ " 1 W 1 -^ i ( 7 ). From the definition of B^(R+, 7 ), we see 

that 

^B m (l+,7)«^/C m (7). 

Suppose now the system (11) has an exponential dichotomy in E m, T^(s) for 
£ G [C -1 ? 1 ] an d Res > 7. An equivalent norm for £"(7) is 

/oo 
|0(s)| m 1 rfw] 1/2 , S = <7 + KJ. 
-00 1 j 

Based on this, one can show that (11) also has an exponential dichotomy in JC m ("/) 
for £ G C}- The definition for exponential dichotomies in a Banach space like 

JC rn (-f) is standard, and can be found in [9]. Using the definitions of the function 
spaces and exponential dichotomies, it is not hard to show the following. (See [13] 
for the case m = 1.) 

Lemma 1. For any <j> G B™ (R + , 7) , consider 

u 2 = K 1 £- i (T i (z,e-,s)p:(e,s)4>(s)), t<?. 

IfsuPteiti-WfiKkKton-ilMOl < ™, then Uj G H 2 ™' 1 ^- 1 x R+, 7 ), 3 = 
1,2 and is a solution to (10) with Uj(£,0) = 0. Moreover 

\Uj\ H 2 m ,l^) < C\<j>\ B ™(R+ n )- 

4. A sequence of functions p G F l , i G Z, where i* 1 * is a Banach space, will be 
denoted by {J 1 }. Define the norm |{/ 4 }|_Fi = sup^l/ 1 !^ »}. 
Consider a linear system, 

2m-l 

(12) uj. + (-l) m D| m « < - ^ A)(0Dlu l =h i (£,T), in£\ 

i=i 

with jump conditions (5), initial conditions (6) and compatibilities (7). Assume 
that the coefficients Aj(£) are extended to £ G R by constants outside fi*. After 
the extension, assume that the associated homogeneous dual system (11) has ex- 
ponential dichotomies in E m - M (s) 1 v = 0, l/4m, for £ G R and Re(s) > 7. Assume 
that g l G i 2 (S 4 , 7 ), 5 i G J B m (R+, 7 ) and uj G H m (Q i ). In this section S 4 = fi* x / 
with / = R + . Assume 

I Wis- (ir+, 7 ) + |{^}|l2(s% 7 ) + |{uo}lff m (0') < 00, 

We look for solutions u l G fl' 2m ' 1 (E i ), i G Z. 

By a standard method, we can continuously extend h l and w to £ G R so 

that 

I^Il 2 (rxr+) < C\h l \ L 2(Y,i), |woleH"*(R) ^ C| u ol-H" m (o*)- 
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We first solve (12) in the domain R x R + , with an initial condition u l but 
no jump conditions. Denote the solution by u l . From the existence of exponential 
dichotomy in E m '°(s), for Rc(s) > 7, we can prove that (12) defines a sectorial 
operator A 1 in L 2 (R). Moreover, 

\\-AT L l m <j^y ReA > 7. 

The proof of the case m = 1 can be found in [13]. The general case can be proved 
similarly. It is than easy to see that u l can be solved by the analytic semigroup 
e A%T and the variation of constant formula. We have 

|ri , |//2m,l( RxR + i7 ) < C(|ft, J | L 2( RxR+ 7 ) + |Uo|jJm). 

Let U i = {u\D(:u\--- ,L>| m - 1 u i ) r , and S l = •) - U i+1 (C, ■). Then 

|* < |Bm(R+, 7 ) < C(|{^}| L 2 (S%7) + |{lij}|ffm(n«))- 

Let the solution to (12), (5)-(7) be u l = u l + v\ The function v' 1 satisfies (12) 
with ti = and i/(£,0) = 0. Let rf = 8 l - 5\ Then rf G B™(R+,j). Let 
V 1 = (v\D 5 v\ ■■■ , D| m ~ V) r . Then the dual systems for V i are 

/ \ 

(13) V£ = J{s)V l + {-l) rn \ ••• \V\ 

\4(o 4(0 ••• 
(w) v i {C,-)-v i+ \c,-) = ff- 

We want to solve (13) and (14) with ff G £"(7). 

I^k-W < C{\5% m[VL+n) + \{V}\ L ^ n) + \{ul}\ Hm ). 

For two subspaces M © N = R 2m ™, denote by P(M, N) the projection with 
the range and kernel being M and N respectively. Assume that at each , we have 

TlP l u (f , s) © ftP] +1 (f , a) = R 2m ™ • 

Here 7£ stands for the range of an operator, and P^ and P l s +1 are projections 
associated to the exponential dichotomies in tt l and respectively. Assume 

that the norms of the projections associated with the above splitting arc uniformly 
bounded with respect to i G Z, Res > 7 in the E m,v (s) norm. Notice that the 
assumption is valid if we choose 7 > to be sufficiently large, due to the Roughness 
of Exponential Dichotomies again. 

We now solve (13), (14) by an iteration method that is used to prove the 
Temporal Shadowing Lemma, [11]. As a first approximation, let 

s) = P{np:{c, s), tzpi +1 {c, s))f(s), 

m 1 - 1 ,*) = -pyip;(?-\8),nptHe- 1 ,s))v i -\s), 



LOCAL AND GLOBAL EXISTENCE OF MULTIPLE WAVES 



7 



From Lemma 1, m^- 1 ^ (£, s)) £ H 2 ™' 1 ^ x M+, 7 ) and is a solution for (12), 
with h 1 = 0. However, at £ l , the jump is not exactly ff . Let •) — </>' +1 (£', ■) = 

ff —f\\. Then 

= r m (f,f +1 , , s) - t\c ,C -\ -\ s). 

Here d = inf{£ l+1 — The above process can be repeated with {fj 1 } replaced by 
{r)\}, and the second approximation denoted by {<p\}. By iteration, we can have 
a sequence {f))}, j > 1 and a sequence of approximations {</>}}• Suppose now the 
constant C\ = Ce~ Q<i < 1, then the convergent series 

oo 

is the desired solution to (13) and (14). 7r 1 £ _1 (<i> 1 ) is an exact solution for v*. The 
uniqueness of {V} can be proved by the exponential dichotomy argument and will 
be skipped here. Observe that by the Paley- Wiener Theorem, 

(15) Tf{ T ) = for r < At, i £ Z, => v l (£_, r) = for r < At, i £ Z. 

This fact will be used in the next section. 

Finally, the solution to the system (12), (5)-(7), u % — u % + v l , satisfies 

(16) \u l \ H 2 m ,i ( ^ n) < C(\{5 l }\ Bm(WL+n) + \{h l }\ L 2 ( ^ n) + \{ui}\ Hm{Qi) ). 

5. The nonlinear system (4)~(7) can be solved by using the result of §4 on the linear 
system and a contraction mapping in a suitable Banach space. The following Local 
Spatial Shadowing Lemma is the main result of Part I. 

Theorem 2. Assume that f is C 00 with bounded derivatives with respect to all 
the variables, {w 1 } is a formal approximation with |{tu*}|j/2m.i(E<) < oo. Let L = 
[0, At] . Assume that e > is small so that a near identity change of coordinates 
can be made in [0, eAr] as in §L Let d = infi{£* +1 — £ 1 } > 0. Then there exist 
f3o, eo > and a positive linear function fJ,*((3), < (3 < [3$. I/O < e < e , and 

IKIIiw + \{g l }\^m + m\ B ^i) < /**(/?), 

then there is a unique solution {u 1 } to (4)-(7), with \{u l }\ H 2 m ,n^ < [3. Moreover, 

|K}| W (E*) < C(\{ui}\ Hm{Qi) + \{g*}\ L 2 m + \m\ B -(i))- 

Proof. Let h % £ i 2 (S 4 ), i £ Z. Since 7 > 0, it is easy to extend the domains of 5 l 
and h % to t e M+ , so that 

|{ft*}|l,2(n*xlR+,7) + |{**}|b"'(R+,7) < C(|{^}|i 2 (Si) + 

Consider the associated linear system (12). From the assumptions, it is clear 
that supj ;i fe |^4fe(0l < 00 • Assume that the coefficients have been extended to 
£ £ M by constants, then from §4, there exists M > such that if M > M then 
(12) has an exponential dichotomy in E mv (s), v = 0, ^ for £ e K and s £ Sg(M). 
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This also implies that (12) has an exponential dichotomy in ]C m {"f) for 7 = M, 
of which the exponential decay rate is a = a(l + \/M). By choosing larger M, 
we have C\ = Ce~ ad < 0.5, where C\ is as in §4. The result in §4 concerning the 
system (12), (5)-(7) is now valid. Let the unique solution be denoted by 

{u i }=^({u} ) },{6 i },{h i }). 

Restricting the solution to the finite interval / = [0, At], in the unweighted 
norm, using (16), we have, 

|jV}| H 2 m ,i (5>) < Ce 7AT |{u J }| ff 2 m ,i (s:2 , xR+rr) 

< Ce^{\W }\ Hmm + \W}\ LH&) + \m\ Bm(I) ). 

Let the solution in that finite time interval / be denoted by 

K} = ^(Kun {ti}). 

Consider Q(J3) = {{«<} : u* e H 2 ™' 1 ^), |K}| ff2m ,i < (3}. Let |{uj}| H m (n * 
+ |{^}|b"(j) + \{g l }\ L -m = 11. For {«*} g Q(/3), we have 

|AA>\ < | 5 J | i2 +C(| U f + e Ar|^|) 

<C(/3 2 + eAr/3 + Ai ). 

Consider the mapping 

{<}=^(K}, {S*},{Af\u\g\e)}). 

We have 

|K}| H2m , 1(Ss) < C(/i + /3 2 + eAr/3). 

Let /? be small such that C(3 2 < \(3. Let fj, and eo be small, depending on (3, such 
that C\i < |/3 and CeAr < |. Then jFj maps Q(/3) into itself. One can also verify 
that if (3 is small, then Tj is a contraction mapping. Therefore, there exists (3q > 
such that J 7 / : Q{[3) — > Q(/3) has a unique fixed point {u 1 }. 

Finally, the solution {w 4 } does not depend on the method of extending the 
domain of {S 1 }, {g 1 } to r e M + . This can be verified by using (15). □ 

Remark. In many formal constructions, d — inf{£* +1 — — > 00 as e — * 0. Then 
the condition Ci = Ce~ Qd < 0.5 is valid if e is small. We do not need to choose 
large 7 = M to make a large. 

PART II. Global Existence of Multiple Waves 

6. The multiple wave solutions constructed in Part I exist only for a short time 
t G [0, At], At — eAr. If {u 1 } is not too large, using the output of the previous 
interval as the input of the next time interval, the process can be repeated to 
obtain solutions in [7 At, (j + l)Af], j = 1, 2, • • • recursively. It is shown, in [13], 
that if a formal approximation is defined for t € R + s, under certain conditions, it 
is possible to obtain global solutions for t e R + . In the second part of this paper, 
we summarize the results in [13]. 
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Although the method should work for some higher order parabolic systems, 
such as the phase field equations, [7], to simplify the matter, we will consider a 
second order system, 

(17) eu t = e 2 u xx + f(u,x,e), x e K, t > 0. 

Assuming by the method of matched expansions , we have the formal series 
for the wave fronts, 

o 

and formal series solutions in the £-th regular and singular layers, 

m 

u m (x,t,e) = J2z j uf e (x,t), 


m 

^(£,M) = £e%f(£,i). 
o 

Here "R" and "S" stand for regular and singular (internal) layers, and £ = (x — 
V%e))/e. 

For convenience, assume the the following Periodicity Hypotheses. 

1. f(u,x + x p ,e) = f(u,x,e); 

2. r) e+i *(t,e) = r] e (t,e)+x p ; 

3. u R V +e '\x,t,e) = u m (x-x p ,t,e); 

4. u s ^)(£,i,e) =u se (£,t,e). 

Here x p > and integer i p > are two constants. The periodicity hypotheses 
ensure that all the estimates obtained here are uniform with respect to layer index 
i. They do not play any other rolls and are not necessary. 

Let < (3 < 1 and let the width of the internal layers be 2e /3 ~ 1 . Define 

The family of curves T l = {(x,t) : x = y l {t)} divides the domain into regions 
S l , i <E Z, where S J is between r i_1 and T 1 . A formal approximation can be 
obtained from the matched expansions, 

\ m {x,t,e), i = 2£-l, 

g<( »V(«,e) >t)e)) i = 2t 



w z (x, t, e) 



u 

Here are the assumptions on w l : 
H 1. There exist C, 7 > such that for all small e and i,feZ, 

\w\x,t,e) - u/(x,oo,e)| < Ce^\ (x,t) G R 2 . 

\f] e (t, e) - V e (oo, e)\ + \D tV e (t, e)\ < Ce~^, t e R 
Here w l (x, 00, e) = lim^^ w % {x, t, e), etc. 
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H 2. There exists a > such that in each regular layer E*, i = 21 — 1, 

Rea{f u (w l (x,t, e),x, e)} < -ct 
uniformly for all (x, t) G E*, ieZ and small e > 0. 

H 3. For an approximation w l (£,t, e) in an internal layer, as £ — * ±oo and e — > 0, 
both iu 1 and dw % /d^ approach the corresponding values of w l+1 or u; 1-1 at common 
boundaries. More precisely, if i = 2£, then for any [i > 0, there exist N, eo > 
such that Cq" 1 > N, and for < e < eo, t > 0, 

\W i (Z,t,e)-W i -\y i - 1 {t),t,e)\<ii, for - e^ 1 < £ < -iV, 
|T^ 4 (£,t,e) -^+ 1 (y i (i),t,e)| < H, for e^ 1 > £ > N. 

Here the function VP = (w l ,w^) is expressed in the stretched variable £ — (x — 
V %e))/e. 

Let £ = £(£) be a function of £ such that 

for |£| < e?-\ 
-e?-\ for £ < -e^-\ 
e^\ for £ > e^ 1 . 

For each t > 0, i = 2£, consider the operator A l {t) : L 2 (R) -> L 2 (R), 

^ (t)u = w ce + A?/ (t, e)u£ + /„ {w i (I t,e),r} e (t,e) + e£, e)u. 

H 4. z = 21, t < 0, has a simple eigenvalue A 4 (e) = eA (i) + 0(e 2 ). The 

rest of the spectrum is contained in {RcA < — a}, a as in H2. Moreover, for the 
limiting operator A l (oo), we have, 

A (oo) < A < 0, uniformly for all i = 2t. 

We look for exact solution u that is of the form w l + u % in each £\ The limit 
u l (x, oo, e) describes the correction to w l {x, oo, e) that yields a stationary solution, 
while u % (x, t, e) — u l (x, oo, e), together with w l (x, t, e) — w l (x, oo, e), describes how 
the solution approaches its limit as t — > oo. Therefore, we define the following 
Banach spaces. For a constant 7 < 0, let 

X(fl x R+,7) = {u : u = 1*1+1*2, mi G L 2 (fi), u 2 G L 2 (9. x R+,7)}. 

I U U( 7 ) = l u lU 2 (Sl) + I u 2|l2(s- 2xR + 7 ). 

X 2 ' 1 ^ x R+,7) = {u:u = i*i+i* 2 , 1*1 G H 2 (n), i* 2 G ff 2 ' 1 ^ x R+,7)}. 

\ u \x 2 ^(-y) = l u l|ff 2 (0) + l u 2|ff2,l (slxR+;7 ). 

It can be verified that for u G x R+, 7) or X 2,1 (f2 x R+,7), the decomposition 
i* = 1*1 + 1*2 is unique. 

For T. 1 = {(x,t) : < x < y l {t), t > 0}, we say that u G L 2 (£*,7), 

H 2 ' 1 (E i ,j), X(S*,7) or X 2,1 (E%7), etc., if u is the restriction of a function u G 
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L 2 (M x K+,7), etc., to the domain The norms are defined by 

Ml 2 (£*, 7 ) = m f{|«lL 2 (RxR+,7)}, 
\ll\x 2 ' 1 = m f{|«U 2 ' 1 (RxR+, 7 )}- 

Let 

x k (j) = {u | u = Ul + u 2 , mi e w\ u 2 g 77 fe ( 7 )}, 7 < o. 
X fclXfc2 ( 7 ) = X fel ( 7 ) x X fe2 ( 7 ), fci > 0, fc 2 > 0. 

Similar definitions can be given to spaces L 2 (flxl, 7 ), if 2 ' 1 (fix J, 7 ), 77 2,1 (£n 
7, 7 ), but not X(Q x J, 7 ) or X fc ( 7 ), if 7 is a finite time interval. 

The main result of Part II is the following Global Spatial Shadowing Lemma. 

Theorem 3. Let {w 1 } be a formal approximation of solutions for (17). Assume 
that the Hypotheses HI H4 are satisfied, and the constants 7 and a satisfy —a < 
— 7 < 0. Then there exist positive constants jo, J 2l £o and a negative constant 
7 = 0(e) satisfying the following properties. Assume that {w 1 } is a formal approx- 
imation near {w 1 }, with 

(18) K -w l \ X 2,i { ^ n) < de n , ieZ, 

for some C\ > and ji > 1. Assume that for the approximation {w 1 }, we have 

(19) I <C 2 e>\ j 2 >j . 

Then for < e < 6q, to any locally 77 1 function uq with \uq — t# l (0)|jji(£» n {T=o}) < 
C 2 e 32 , there exists a unique exact solution to (17) that satisfies u exact (x,0) = 
u (x), and 

\Uexact - W l | X 2,l( S%7 )) = 0(e j3 ), i g Z. 

77ere js = min{ji, j 2 — J 2 }, J 2 > is a constant that does not depends on e. All 
the norms in this lemma are expressed by the stretched variables £ = x/e, t = t/e. 

Remark. (1) The constants j and J 2 depend on {w 1 }. However the approxima- 
tion may not satisfy (19). By adding higher order expansions, the new approxima- 
tion {w'} will satisfy (18). Therefore, it has the same constants j and J 2 , which 
are stable with respect to perturbations. Moreover, if the order of expansion is 
sufficiently high, {w 1 } will also satisfy (19). 

(2) The values of jo and J 2 are not important in applications, since j'3 = ji if 
72 is sufficiently large. Computing ji is relatively easy. One only need to compare 
{w 1 } with the next approximation. 

(3) Hypotheses H1-H4 are consequences of hypotheses used to construct 
formal expansions for system (17). Sec §11 for a brief discussion. 

7. To prove the theorem, we use a partition of the time interval t G M + , 
K + = Uplb'Ai, (j + l)At] U [t f , 00), 

where tf = rAt. We assume that At = eAr, where Ar is independent of e. tf is 
large such that e~~ (t = 0(e 2 ), where 7 is as in HI. If e is small the variations of 
boundary curves x — y i_1 (i) and y l (t) are small such that a near identity change 
of coordinates can be made in each £*• = n {t g [jAt, (j + l)Ai]}, j < r — 1 or 
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= n {t e [t/.oo)} to straighten T 4-1 and P. In the sequel, we assume that 
P = {x = y l j} in the jth time interval, where j/j- is independent of t. 
Define 

e=N-(yr 1 +y-)/2]A, inSj., 
L)(e) = (y i j -yi- 1 )/(2e), 
n) = [-L)(e),L)(e)], 

P = [0,At], <j < r- 1, and J r = [0,oo). 

Let r = (f - jAt)/e. Then Ej = fij x P. In each E*-, j < r - 1, we linearize 
around w l at the fixed time r = 0, while in E*, linearize around w l at the fixed 
time r = oo. Denoting the solutions in E* by uj, we have 

(20) ^ = ^ + ^(0^ + 4(04+^(4^,^,6), inE}, 

(21) W j (L i (e),T)-U} +1 (-L i + 1 (e),r)=5 i j (r), at P, 

(22) 4 +1 (£,0) = 4(6, At), 4(6 0) =«j(0. 

Here £i = 6* (0 is a near identity change of coordinates induced by the coordinates 
change that straightens the boundaries in each EJ . V- is the constant wave speed 
in an internal layer, but is slowly £ dependent in a regular layer due to the fact 
that the boundaries are not parallel there. 

The idea of the proof is presented in the next three sections. In §8, we study 
the spectral properties of the linear systems and exponential dichotomies asso- 
ciated to the dual systems of the linear systems. In §9, we study the system in 
the last interval [i/,oo). There we need to obtain asymptotic behavior of the so- 
lution {u l r }. Thus the space JT 2,1 (Xy,7) will be used. The result obtained there 
also serves as an upper bound of the accumulation error in the first r intervals 
[7 At, (j + l)At], < j < r — 1. In §10, we study the system in finite time intervals. 
Since r is large, the estimates obtained in Part I is not accurate enough. More 
precise estimates will be derived which rely on the Hypotheses HI H4, while the 
estimate in Part I does not. 

8. Since nonlinear systems can be solved by contraction mappings, we study a 
linear system first. Drooping the subscript j, we consider, 

(23) < = <4 + ^(04 + A'tOu' + ti(£, r), 

(24) U l (L\e),r) - U i+ \-L i+ \e), r) - 5\r), 

(25) u'&O) =4(0- 

The homogeneous dual system of (23) is 

Here U — (u,u^) T . It is crucial to study the exponential dichotomies for system 

(26) . 
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We first discuss the system in regular layers. Freezing £ = £o and consider (26) 
with constant coefficients. Using H2, we can verify that for Re(s) > —<7q, <to = 
<t/2, the system is hyperbolic, having n eigenvalues with positive (negative) real 
parts. 

Now recall that in regular layers, A l (£) — fu(w t ((y' l ~ 1 + y l )/2 + e£)i (f* 1 + 
y l )/2 + e£, e) depends slowly in £, i.e., dA g^ — 0(e). Also, from the change of 

coordinates, cf. [13] for details, V*(£) also depends slowly on e in regular layers. 
Then a theorem from [3] indicates that (26) has exponential dichotomy for £ e R 
and Re(s) > — uq. 

In [13] it is shown that if (26) has an exponential dichotomy in R 2n for £ e M 
and for every s e {Re(s) > — <r }, then it has an exponential dichotomy in the 
space E 1 ' v (s),v = or 0.25. 

Let the right hand side of (23) define a linear operator A % {t), t = jAt, < 
j < r — 1, or t = oo, j = r. Using the fact that (26) has an exponential dichotomy 
in E 1 ' ,y (s), v = for ^gl and s € {Re(s) > — <jq}, we can show that in regular 
layers, A l (t) is an exponentially stable sectorial operator in L 2 (R), i.e., 

Re{aA l (t)} < (Jo , 

for every t > 0. 

We now discuss (26) in internal layers. Here the spectrum of A l (t) is given by 
H4. Also D t r] e {t, e) = Vj, i = 2l\it = jAt, < j < r - 1, and V* = 0. From H3, 
if fi is small, then for |£| > N, the coefficients A l (^) and V^(£) are close to that of 
A l (±L l (e)) and V i {±L i (e)). But the system with constant coefficients A i (±L i (e)) 
and U l (±L l (e)) is hyperbolic. Therefore, from the Roughness of Exponential Di- 
chotomies, which is also valid in spaces E 1 ' v (s),v = 0, 0.25, (26) has exponential 
dichotomies in E 1 ^(s), v = 0, 0.25, for Re(s) > -er and |£| > N. 

The exponential dichotomies in the region £ < — N and £, > N extend to 
£ G R~ and £ G R+ respectively. Let the projections be P l s +P^ = I. We can show 
that 

(27) ^P„(0~, s) "KP s (0+, s) = K 2 ™, if s 7^ A*(e) and Re(s) > -a . 

If (27) were not true, we could find a solution {/(£, s) for (26) that decays expo- 
nentially as £ — > ±oo. Then s would be an eigenvalue with t/(£, s) as an eigen 
function. This contradicts to the fact A 4 (e) is the only eigenvalue in the region 
Re(s) > —co- 
in [13], it is also shown that the projections defined by the splitting (27) is 
0(1 + s _^( e )| ). Observe that H4 also implies that for the linearization at t = oo, 

(26) has an exponential dichotomy in £ e R in the region Re(s) > where 
A < is as in H4. 

9. We now study (23)-(25) in the time interval [tf, oo). The procedure of solving 
them is similar to that used in Part I. The right hand side of (23) defines a sectorial 
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operator *4 l (oo) in L 2 (R). Extend the domain of u and h % to £ e R. Let 

u l = e Ai ^ T ul + f eA'M^-Qh'iQdC. 
Jo 

As shown in §8, 

ovf(oo) C {Rc(s) < -Ce < 0}, 

both in regular and internal layers. Because of the spectra of A l (oo), with some 
7 = Ce < 0, we can show 

(28) l^l^w < C(e-°- 5 \4\ H1 + e-^I^Uw). 

The solution {u*} for (23)-(25) can be written as u l = u l + v l , with {v 1 } 
satisfying 

(29) 4=^ + ^(04+^(0"*, 

(30) V^(L*(e),r) - ^ +1 (-^ +1 (e), r) - iftr), 

(31) «*(£,<)) = 0. 

Here F = (v, v^) T , 77* includes <P and corrections from the jumps of {u 1 } at bound- 
aries. We have 

177^0^x0^(7) < C(|(5 4 | x o.75xo.25 (7) + e-°- 5 \{u }\ H i +e- 1 - 5 |{^}| i 2 (Ei;7) ). 

Since for the dual system of (29)-(31), the exponential dichotomies exists in the 
region Re(s) > 7, 7 < 0, using the Laplace transform and the iterative scheme 
similar to that used in Part I, we can show that system (29)-(31) has a unique 
solution {v 1 }, with 

KU 2 .!(7) - C\W}\X°™X0 .25 ( 7 ). 

Therefore, in the last interval, the solution {u 1 } satisfies 

\u l \x^ {l) < C(|{<5 4 }| x o. 7 5xo.25 (7) + e-°- 5 \{u }\ Hl + e- 1 - B |{/i < }| L2(E4i7) ). 

Using the above result, the nonlinear system can be solved by a contraction 
mapping in a ball of radius e ri , 77 > 1.5, in X 2 ' 1 ^). The following result is proved 
in [13]. 

Theorem 4. Let \g*\ x{y) = (e^+ 1 - 5 ), |uj| H , = o(e^+ - 5 ), 1^^0.75x0.25^ = 
o(e ri ), where 77 > 0.5. Lei be e dependent such that e~ 7t/ < C$e 2 , where 7 
is as in HI and C is independent of e. Then there exists e > such that for 
< e < e ; the nonlinear system (20)- (22) has a unique solution {u 1 } satisfying 
the following estimate, 

(32) |K}U2. 1(7) < C{|{5 i }| H o.7.xo.» (7) + e-°- 5 |{4}| Hl + e- 1 - 5 !^}^^)}. 

10. We study system (20)-(22) in the finite intervals, [7 At, (j + l)Ai], < j < r-1. 
The existence of solutions in a finite interval has been discussed in Theorem 2. 
However, to guarantee the existence of solutions up to the last interval [tf, 00), a 
stringent restriction on the accumulation errors in these finite intervals must be 
met. Since the critical eigenvalue A'(e) may not be negative, u l - may grow as j 
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increases. From Theorem 2, {uj} is determined by {m*(0)}, {5*} and {Sj}. Among 
them only {uj (0)}, carries the information from the previous j intervals. Therefore, 
it is crucial to control |u*-(Ar)| in terms of |{uj(0)}|. 

Let Qq J and Q % j be the projections corresponding to the spectral set in 
A'(e) and {Re(s) < — 00} in the ith internal layers and the jth time interval. Let 
Qo 3 Uj = C(j(f)j where <f>j is the eigenvalue corresponding to the eigenvalue A l (e). 
Then 

\a)\ + \Q?u)\ H i, 

is an equivalent norm for The following estimate has been proved in [13], 

|aj(Ar)| + \Q^u)(Ar)\ m < (1 + Ce) sup(|a}(0)| + \Q^u)(0)\m) 
(33) 1 

+ C|{^}|// - 75 x°- 25 (/j) + C\{g l j}\ L 2. 

Since the spectral projections for the j and (j + l)th interval differ by 0(e), 
K+iWI + IQs'-'+X+iWIhi is also bounded by the right hand side of (33). There- 
fore, one can show that 

\a)(0)\ + \Q^ui(0)\ m < (l + Cey'sup(|aj(0)| + |Qi>°t4| ff i) 

i 

+ ( 1 + Ce ) J SU p(|{^}| g0 .75x0.25 + \{gi}\ L ')- 

e k<j 

Using the fact (l+Ce) 1/(Ce) < e and j < r w log(i)/e, we see that (l+Cep < 
e~ B for some B > 0. Therefore, if all the terms in the right hand side are bounded 
by e M for some large M > 0, have the solution for all the finite intervals, and 
the initial condition for the infinite interval is small enough so that Theorem 4 can 
be used there. Here is the main theorem in this section. 

Theorem 5. There exists At > such that if r = [log(— ^)/e7Ar] + 1), then 
we have the following concerning the solutions of (20) -(22) in P , < j < r — 1. 
There exist M\ and M2 > 0, such that if 

\{4(0)}\ m + \{5*}\ H a, 5 «a, 5 + \{g*}\ L2 = 0(e M >), 

uniformly for < j < r — 1, then (20)-(22), < j < r — 1 has a unique solution 
for < e < €0, where e > is a small constant. The solution in each P satisfies 

|aj(0)| + \Q^u)(0)\ m < Ce- M2 [su P (K(0)| + |Qj°«j| ffl ) 

i 

+ SUP(|{^}| H 0.75X0.2 S + |{5fe}U 2 ]" 

k<j 

In particular, if Mi is sufficiently large, in the rth (infinite) interval \{u l r (0)}\ H i = 
o(e ri+0 ' 5 ), r\ > 1.5. Therefore, Theorem 4 applies in the last infinite interval also. 

To describe the idea of the proof of (33), it suffices to consider the linear sys- 
tem (23)-(25). For simplicity, we consider the Oth interval and drop the subscript 
j. The procedure of solving such system is similar to that used before. Define u l as 
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in Part I, and let u'(t) = a^r)^' + Q*'°m*(t). In internal layers, A* = 0(e). Also 
in the stable space, the semigroup is exponentially stable. We can show, 

\a l (Ar)\ < e ( 1+Ce ) AT (K(0)| + VAr^V), 
\Q l fu\Ar)\ m < C(e~^\Ql- a u\0)\ m + \a\0)\), 

for some <r < 0. Observe in the above, u l (0) = u l (0). Thus a 4 (0) = a*(0). 

Let At > be large so that C'e- aAr < 1/2. We then find d such that 

e (l+Ce)Ar < x + (j ie _ We haye 

(34) |^(At)| + \Q l fu\Q)\ m < (1 + de)K(0)| + l/2|QjV(0)| ff i + C|^| L2 . 

The above is also valid in regular layers with a 1 = 0. 

The solution for (23)-(25) again is written as u l — u l + v l , where {v 1 } sat- 
isfies (29)-(31). As in Part I, system (29)-(31) is solved by iterations, and {v 1 } 
is bounded by {rf} in suitable function spaces. An complication arises here since 
{rf} also include corrections of the jumps of {u 1 } at common boundaries. The 
latter in term is bounded by |{a*(0)}| and |{Q*'°u 4 (0)}|, of which the coefficients 
have to be carefully controlled. 

It can be shown that {rf} depends weakly on a*(0), since \a t (j) l \ < Ce _cr ^ 
for some a > 0, and L z (e) — > oo as e — > 0. However, it is a little tricky to control 
the dependence of on |{Q*'°m*(0)}|. The key here is to use the exponential 

decay of e A T in the stable subspace (which is the whole space in regular layers). 
To do so, weighted norms have been used, [13], 

|C/|fl"0.75X0.25(/J iO .) = \e <TT C/|/J0.75x0.2S(-/j), U < 0, 

\ u \h 2 - 1 {WxIi ,a) = \ e ^' 7Tu \H 2 . 1 (n i xD)i a < °- 

Observe here the interval P = [0, At] is finite. It is shown that 

|?7 J |#0.75x0. 2 5(/ 3jO .) 

+ sup{e|a l (0)| + C|<%V(0)| + Ce^ AT \h l \ L 2}, 

i 

for some a < 0, independent of e. Also it is shown that the solution of the system 
(29)-(31) satisfies 

|w*|H2,l(Qi x 7J,(7) < C|{77*}|i/0.75X0.25(/3 iO .). 

Let QjV = a>\ 

W(At)\ + |£#V(Ar)| ff i < C|^(AT)| ffl 

< Ce aAT \v l \ H 2,i^ XIJ ^ 

< C|{(5 I }|jj0.75x0.25(-/3 j 

+ sup{Cie <rAT |Qj J '« i (0)| jr i + Cee aAT \a*(0)\ + C\ti\ L2 }. 
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Let At be large so that Cie crAT \Q 1 ^ u l (0)|#i < 0.5. Combine the above with 
(34), we have an estimate for the linear system, 

K(Ar)| + \Q^u\Ar)\ m < (1 + Ce)sup(K(0)| + |QjV(0)| ff i) 

(35) 

+ C|{^}| H 0.75X0. 2 5 (/ , ) +C|{/ l l }| L 2. 

The estimate of the nonlinear system (33) can be obtained from (35) by replacing 
h l with Af l , and using the fact that if |u l | ff 2,i < e, 

I^V < W\ L -+Ce\rf\ H ^. 

11. Constructions of multiple waves in singularly perturbed equations have been 
discussed in many papers, see for example [5]. Expansions for system of reaction- 
diffusion equations to any order of e can be found in [12]. The following hypotheses 
are used in that paper. 
There is a partition 

R = UZ- 00 [x e ,x e+1 ] 

that is periodic with respect to £, compatible with the period of /. A C°° function 
p l {x) is defined on [a; 1-1 , a;*] with fo{p % (x), x) = 0. Also assume that p e+ip (x+x p ) — 
p^(x), x G [x 1 ^ 1 , x^\. 

H* 1 Rca{f 0u (p l (x),x)} < for x e [a;*" 1 ,^], i G Z. 

— ■ x 1 

Using a stretched variable £ = , we assume that the 0th order expan- 

e 

sion of Eq. (17), 

(36) u^ + fo(u,x l ) -0. 

has a heteroclinic solution <?*(£) connecting p l (x l ) to p 4+1 (x 4 ). Assume 
H* 2 The linear homogeneous equation 

0« + /o u (g i (O.^ = o, 

has a unique bounded solution <7|(£), up to constant multiples. 
From H*2, we infer that the adjoint equation 

+ /on (q\{0,x i )1> = o 

has a unique bounded solution V>i(0 U P to constant multiples. 

The following assumption implies that the heteroclinic solution breaks as x 
moves away from x 1 . 

H* 3 tfiOfoxtfi&x'W ? 0, 

When a; is in a neighborhood of x 4 , we look for traveling wave solutions with 
wave speed V l {x). It is also of interest to find out conditions to ensure that the 
wave front moves towards x % . Consider A 1 : L 2 (R) — ► L 2 (M), defined as 

A'u = u ss + / „(<f (0,*>, D(A l ) = H 2 (R). 
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H* 4 A = is a simple eigenvalue for A 1 , ieZ. There exists ao > such that 
(j{A l ) C {A = 0} U {ReA < -cr }- 
Since A = is simple, we have 

/oo 
-oo 

h* 51 = [j^ i&no^mr 1 /-°°oo ^(o/te^'K). > o, » g z. 

Under Hypotheses H*l H*5, it is shown, in [12], that formal expansions of 
wave front positions {rf(t, e)} and matched expansions of solutions {w l (x,t, e)} 
can be constructed. These hypotheses are general in the sense that H*2 and H*3 
are generic assumptions. Hypothesis H*l is equivalent to the stability of p % (x) in 
a regular layer, as an equilibrium of an ODE, obtained from the 0th expansion of 
(17), i.e., by setting eu xx = 0. Hypothesis H*4 is equivalent to the stability of the 
traveling wave solution q l in an internal layer, in the sense of Evans [4] . Hypothesis 
H*5 implies that the wave front is moving towards x l , if x is near x % . In fact, if is 
shown in [12] 

^--'•<—>- 

Similar conditions like H*5 was used by Fife [5]. Without such conditions, the 
wave front may move away from the stationary position and the formal solutions 
may not exist globally in time. 

It is shown in [13] that Hypotheses H*l H*5 imply Hypotheses HI H4 of 
this paper. Only H4 needs some explanation. When {w 1 } is a 0th order approx- 
imation, A = is always an eigenvalue. When higher order terms are added, the 
eigenvalue becomes O(e). In [13], a precise formula expressing A*(e) as a function 
of the wave speed and the wave front position was proved. In particular, when 
t = oo, it yields 

d\ l _ dV i 
de dx 

Here Vx) is the wave speed as a function of the wave position x. It is now clear 
that H*5 implies H4. 

Many authors have found that the signs of A l (e) or the derivatives of the 
wave speed determine the stability of the multiple waves, [6]. The precise relation 
between the two quantities in general systems is interesting in its own right right. 

In conclusion, the formal approximations obtained in [12] satisfies our hy- 
potheses in Part II. Therefore, the Global Spatial Shadowing Lemma can be used 
to ensure the existence of a global solution near the formal approximation. 
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